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Large neutron star masses as the recently measured 1.97 ± 0.04 M© 
for PSR J1614-2230 provide a valuable lower limit on the stiffness of the 
equation of state of dense, nuclear and quark matter. Complementary, the 
analysis of the elliptic flow in heavy ion collisions suggests an upper limit 
on the EoS stiffness. Wc illustrate how this dichotomy permits to constrain 
parameters of effective EoS models which otherwise could not be derived 
unambiguously from first principles. 
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1. Introduction 

Neutron stars (NS) which outrange the domain of 'typical' NS masses 
(very roughly between 1.2 and 1.6 Mq, see Fig. 1 in [1] for an overview) 
have been considered to provide a serious constraint on the stiffness of the 
equation of state (EoS) [2] but treated rather cautiously due to the fact 
that in the rare cases of very massive observed NS either the accuracy of 
the measurement has been worryingly low or /and the measurement itself 
raised doubts within the community [3]. It seems, that this situation has 
changed after a mass of M=1.97 ± 0.04 Mq has been reported for PSR 
J1614-2230 with an unprecedented accuracy in this high mass regime [1] 
and without perceivable objections from expert groups involved in this field. 
This observation of a two solar mass NS provides a very direct constraint 
on the minimum stiffness of the EoS of cold and dense matter and there- 
fore promises new insights regarding our understanding in particular of the 
nature of particle interactions at finite densities. Implications have been 
discussed shortly after the measurement became public [T]. 
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As an example (which deserves more attention than can be granted in 
this article) for the importance of this result it is worth to mention the 
problem of strangeness for the hadronic EoS. Even though we are rather 
interested in the possible existence of a quark matter (QM) core in NS one 
can not exclude the possibility that before the according critical density 
can be reached the mass-density threshold for hyperons is passed and the 
EoS has to soften due to the appearance of these additional degrees of free- 
dom. Due to this softening the maximum mass of a corresponding NS can 
decrease drastically in comparison to the plain EoS with only neutron and 
proton degrees of freedom and is therefore not unlikely to contradict the 
high NS mass constraint. The problem can be avoided by accounting for 
repulsive vector interaction terms which are a specific feature of relativistic 
approaches to EoS. This has been shown, e.g., in a generalized nonlinear 
Walecka-type model [SJ [6] and a quark-meson-coupling model [3 [8]. Simi- 
larly, it turns out that vector interaction terms are of crucial importance if 
aiming at the description of high mass NS with QM cores [9]. It is an in- 
teresting finding that merely by the observation of a two solar mass NS the 
vector interaction can be identified as an inevitable channel if approaching 
a microscopical description of both, dense hyperon and QM. 

Besides a two solar mass constraint limiting the 'softness' of the EoS we 
suggested several other constraints from NS observations and the analysis 
of heavy ion collisions which have to be fulfilled simultaneously by a viable 
state-of-the-art EoS |2]- In detail, the scheme suggests that a viable EoS 
should 

• reproduce the most massive observed neutron star, 

• avoid the direct URCA (DU) cooling problem, 

• result in neutron stars within the predicted mass-radius domains of 
4U 0614-1-09 (deduced from quasiperiodic brightness oscillations) and 
RX J1856-3754 (deduced from the objects thermal emission), 

• explain the gravitational mass and total baryon number of pulsar PSR 
J0737-3039(B) with at most 1% deviation from the baryon number 
predicted for this particular object, and 

• not contradict flow and kaon production data of heavy-ion collisions. 

At the time of publication of Ref. [2], the most massive NS has been PSR 
J0751-I-1807 with M ~ 2.1 Mq, a result which later has been withdrawn 
[lOj . We reasoned, that in particular the flow constraint is an extremely use- 
ful constraint. In contrast to the two solar mass constraint it limits the EoS 
such, that it cannot exceed the upper limits on the pressure (as a function 
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of density) in symmetric matter as obtained from the analysis of the ehiptic 
flow of iso-symmetric matter in HIC [11]. This gives two complementing 
phenomenological constraints on the EoS stiffness, an upper (flow) and a 
lower (two solar mass constraint) limit, which both considered together sig- 
nificantly reduce the possible shape of the pressure density relation, viz., the 
EoS at super saturation densities. In previous work we took advantage of 
this insight and adjusted otherwise not well determined coupling constants 
of a Nambu-Jona-Lasinio (NJL) type model. Details about this models 
framework is found in the original work formulated before the invention of 
the quark model of elementary particles \\.2\ [T3] and a number of review 
articles [HI [151 [13 [13 [IB] applying it to elucidate the role of dynamical 
chiral symmetry breaking for hadron structure and its restoration in dense 
quark matter. Our resulting hybrid, nuclear-quark matter EoS, is in agree- 
ment with both of these before mentioned constraints and shows a better 
overall performance than the originally underlying nuclear EoS based on the 
Dirac-Bruckner Hartree-Fock (DBHF) approach [9], which will be shortly 
summarized later in this paper. At this time we had to scan a small pa- 
rameter range only in order to obtain this result. However, the question 
how well the QM model EoS is constrained, viz. which other sets of cou- 
pling constants would reproduce a phenomenologically sound EoS, has been 
left open in [21 • The reporting of PSR J1614's high mass has triggered 
this kind of more systematic studies for Bag-like QM EoS within a purely 
phenomenological model for the QM equation of state consisting simply of 
a power series expansion in the quark chemical potential /x including 4*^* 
and 2"'^ order terms [20] . The 4*'*-order term is thought to mimic the 
influence of strong interactions on the ideal gas expression for the quark 
pressure, while the 2"*^ order term is a measure for the competing effects of 
a finite strange quark mass and a possible diquark condensate in deconfined 
matter |21j . It is worth to notice that a ^^-term is not necessarily the lead- 
ing order term in a /u-expansion of the pressure, as has been shown for a 
simple, semi-analytic model based on Dyson-Schwinger techniques at finite 
densities [22]. We would like to mention a very promising, recent develop- 
ment in modeling quark matter in the nonperturbative, low-energy domain 
of QCD which is relevant for the QCD phase diagram and compact star 
phenomenology. This concerns nonlocal, separable interaction models with 
either covariant or instantaneous formfactors [23[ [2^ [25] and their applica- 
tion to compact stars |26[ [27] generalizing the local current-current coupling 
of the NJL model. In particular the rank-2 separable models which allow 
a simultaneous description of the dynamical quark mass function nn{p) and 
the wave function renormalization Z[p) of the quark propagator in accor- 
dance with lattice QCD data [281 (221 [Ml [M] shall allow to greatly reduce 
ambiguities of the parametrization of quark models discussed in these Pro- 
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ceedings. For the time being, however, the nonlocal models need to be 
further improved to address diquark condensation and asymmetric matter 
before they can be used to study compact star constraints. 

In this paper we will illustrate in a systematic, 'whole-range' scan of an 
NJL model EoS for QM how strongly a two-solar NS mass measurement 
constrains the strength of available (vector and diquark) coupling strength 
parameters and what implications can be derived for the investigation of 
HIC, viz. the EoS of iso-spin symmetric matter. Section [2] gives a summary 
of the applied QM EoS, section [3] discusses the result of our analysis, section 
Hlwill provide conclusions and a brief outlook concerning further interesting 
questions. 

2. Hybrid Matter EoS based on a NJL Model 

In order to obtain a QM equation of state we employ a three-flavor color 
superconducting NJL model with selfconsistently determined quark masses 
and diquark gaps [32l [33l [M] . In addition to a typical scalar interaction 
term, we account for a repulsive vector interaction term which stiffens the 
EoS with increasing interaction strength and results in sufficiently high NS 
masses (details are found in [9] and references therein). The attractive scalar 
diquark channels are responsible for the formation of diquark condensates 
and color superconducting phases in the system. As we discussed before 
it moreover lowers the transition density to a nuclear matter EoS with in- 
creasing coupling strength [9]. The effective Lagrangian can be split into a 
free particle and an interaction part. The free part reads as 

Ckin = qi-ilf^df, + m + fi)q, (1) 

where m = diag(m„, m^, m<() is the diagonal current quark mass matrix 
and p, the corresponding quark chemical potential matrix. The effective 
interaction is written as 

8 

Ant = GsrjD ^ {qil5TahCq'^){q^ Cij5TaXaq)+Gs ^ [{qTaq)"^ + Vviqijoq)'^] ■ 

a,b=2.5J a=0 

. (2) 

Ta and Xa are Gell-Mann matrices in flavor and color space respectively and C 
the charge conjugation matrix. Here we have omitted interaction channels which 
do not contribute to the thermodynamics at mcanfield level, like the pseudosclar 
isovector channel which would be required to make the chiral symmetry of this 
interaction model manifest. The parameter Gs defines the scalar coupling strength 
and can be determined from meson properties in the vacuum. For this study 
we apply the parameters obtained for a sharp cut-off regularization scheme from 
Table III of Rcf. [35] labeled with 'oo' in front of the corresponding table row 



KlahnCS-HIC-QCD printed on December 1, 2011 



5 



Q. The parameters rjjjirjv) are defined as the ratio of the diquark(vector) and 
scalar couphng. Since these two parameters are not fixed by vacuum properties of 
mesons or hadrons we treat them as free model parameters. In order to investigate 
thermodynamical properties of the system we use the partition function in path 
integral representation, 



Z{T, {i) = j VqDq exp Ij^^dTj (fx[q{i$ - m + lLf)q + C^^t] \ ■ (3) 



Bosonic meson field degrees of freedom can easily be introduced by applying cor- 
responding Hubbard-Stratonovich transformations in all interaction channels. For 
the sake of simplicity all of them are treated on the meanficld level, viz., mesonic 
fluctuations and higher correlations are neglected. Minimizing the thermodynami- 
cal potential = — T In Z with respect to the meanfields then defines the pressure 
p = — f2 of the equilibrated system. As a result of the minimization procedure 
one obtains a set of coupled gap equations which has to be solved selfconsistently. 
Finally, the thermodynamical potential reads as 



n(T, n) 



SGs 8G\/ 
d^p 



/ E + ^Tln (l + e-^"/^)] no. (4) 



(27r)3 ^ 

^ ' n=l 

Qi denotes the contribution of electrons and muons to the thermodynamic po- 
tential, flo is the contribution to be subtracted in order to obtain zero vacuum 
pressure. The 0/ (with / = u, d, s) are the chiral condensates corresponding to the 
three quark flavors which are obtained from the previously described minimization 
of the scalar meanfleld terms. Accordingly, on has to account for the vector mean- 
flclds Lj f and pairing gaps Ayy . Under neutron star conditions one additionally 
imposes electric charge neutrality and /3-equilibrium conditions. In this paper we 
disregard the Kobayashi-Maskawa-'t Hooft term in the interaction Lagrangian, see 
Ref. |33| for a motivation. Further references and possible consequences arising 
from the inclusion of this term are discussed in [321 ISZ] • 

With Eq. (HI we have a QM EoS available which is still, by construction, phe- 
nomenological but allows to interpret the influence of certain interaction channels 
on the EoS more specifically than generalized power expansions of ideal gas ex- 
pressions. The next step in order to describe NS phenomenology is to account for 
confined nuclear matter at lower densities. Without a unified in-mcdium-approach 
for the description of nuclear matter in terms of QM degrees of freedom at hand 
we, as everybody else, fall back to a two-phase description, joining independently 
obtained nuclear and QM EoS by performing a phase transition construction. As 
common as this procedure is we feel a few comments have to be made about it in 
order to avoid a misleading interpretation of this work. Our understanding of the 



^ This parametrization scheme has been implemented in an online tool developed by 
F. Sandin which also corrects for a mistake in the kaon mass formula employed in 



[55] , see http://3fcs.pendicular.net/psolver 
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nuclear EoS at supersaturation densities is by no means more profound than our 
knowledge about dense QM. The set of constraints we apply in order to pin down 
the QM EoS we originally bundled up in order to constrain the nuclear EoS. There 
are no profound additional insights into the shape of the 'true' nuclear matter EoS 
since the publication of [5] even though progress has been made in order to gain 
deeper insights into the physics of finite density nuclear systems, e.g., within the 
framework of chiral effective field theory [35] and directly from QCD via lattice 
simulations [39]. While these studies are progressing and may become applicable 
to high density nuclear systems in compact stars and heavy-ion collisions in future, 
we discuss these systems for the time being from the point of view of parametric 
approaches to the high density EoS. Such cold dense EoS studies may be guided 
by observations of compact stars, see [U UOl EI]- From this parametric point of 
view the variation of model parameters of a QM EoS in a strict sense would always 
require to perform a similar variation of an independently obtained NM EoS. All 
results of this (or any similar) analysis, are likely to change significantly if the NM 
EoS is exchanged by a model with significantly different high density behavior. As 
an example, softening the NM EoS to a degree where the model is not any longer 
in agreement with the two-solar mass constraint requires to soften the QM EoS as 
well in order to obtain a thermodynamically sound phase transition. This makes 
it hard, if not impossible, to obtain a hybrid EoS which would describe a two solar 
mass NS. Still the QM EoS by itself is not necessarily in contradiction with the 
existence of high neutron star masses. 

For our analysis we avoid the problem of a nuclear parameter scan by applying 
the Dirac-Brueckner Hartree-Fock (DBHF) EoS which has proven to perform rea- 
sonably well for describing nuclear matter saturation properties and kaon data |42j 
as well as NS properties [2] even though it tends to behave too stiff above densities 
of about 3.5 times saturation density. On the other side, this stiffness occurs in a 
region where QM degrees of freedom are not unlikely to be the only ones which are 
relevant. Amongst other reasons we prefer the DBHF EoS because it is based on 
a relativistic and microscopical description of many-particle interactions. It starts 
from a given free nucleon-nucleon interaction (the relativistic Bonn A potential) 
fitted to nucleon-nucleon scattering data and deuteron properties. In ab initio 
calculations based on many-body techniques one then derives the nuclear energy 
functional from first principles, i.e., treating short-range and many-body correla- 
tions explicitly. In the relativistic DBHF approach the nucleon inside the medium 
is dressed by the self-energy based on a T-matrix. The in-medium T-matrix as 
obtained from the Bethe-Salpeter equation plays the role of an effective two-body 
interaction which contains all short-range and many-body correlations in the lad- 
der approximation. As we have shown in the context of hybrid EoS the rather stiff 
behavior at high densities is not necessarily relevant if the phase transition to QM 
occurs at low enough densities of about three to four times saturation density [5]. 

3. Results 

3.1. Hybrid Neutron Stars 

Unlike to previous work, where we applied both, the flow and two-solar mass 
constraint, simultaneously in order to obtain a phenomenologically sound hybrid 
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EoS [S] we start this analysis from calculating neutron star configurations for a 
wide range of vector and diquark couplings [rjy £ [0.0,0.7], rjo G [0.8,1.15]). To 
further extend the previous study we apply a different parameterization for the 
scalar coupling strength, as well, applying the parameters from the row labeled 
'cx)' in Table I of Ref.[3S]. As we will show, this choice significantly affects the 
outcome of the present study. 

Before we discuss the overall result of a full variation of the two free parameters 
rjv and rjD we keep one of them fixed at a reasonable value and vary the other. We 
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Fig. 1. At constant vector coupling (here rjv = 0.3) an increase of rju lowers both, 
the critical density and the maximum NS mass. 



consider any choice to be reasonable which eventually describes a two solar mass 
NS if one varies the one remaining free parameter. In Fig. [T] the vector coupling 
is kept constant at a value of riv = 0.3 while the diquark coupling is varied in the 
range r]D = 0.80 . . . 1.10. At 'low' values of rin, here up to r]D = 1-00 we find the 
required massive NS configurations. As one observes, an increase of rju does not 
only result in a decrease of the maximum NS mass but lowers the critical density for 
the phase transition, too. In other words, an increase of r]D increases the content 
of QM in massive NS and lowers the maximum NS mass at the same time. On the 
other hand, keeping the diquark coupling 77/3 at a constant value and increasing the 
vector coupling will increase both, the maximum NS mass and the critical density. 
This is illustrated in Fig. [51 Again, there is a critical value of rjv corresponding 
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Fig. 2. At constant diquark coupling (here, rjjj = 1.0) an increase oi rjv increases 
the maximum NS mass and the critical density. 



to a minimal stiffness of the EoS which has to be exceeded in order to obtain NS 
configurations with a sufficiently high maximum mass. In the illustrated example 
for 77/3 = 1.0 this holds for values larger than a critical coupling slightly above 
VD = 0.2. 

With this understanding of the influence oi rjo and rjv on both, the maximum 
NS mass and the critical density for the onset of the phase transition we now 
perform a variation over the full available parameter space region of rjD-rjv in 
which we can obtain stable hybrid NS configurations. The result of this study 
is summarized in Fig. [3J which we will now discuss in detail. From the obtained 
mass-radius and mass-density relations of each of the differently parameterized 
hybrid EoS we extracted two numbers only, the maximum possible NS mass and 
the NS mass at which the central density becomes large enough to generate a QM 
core in the NS. Consequently, we call the latter quantity Monsct- The red band 
in Fig. [3] labeled PSR J1614-2230 corresponds to EoS parameterizations which 
describe a maximum mass exactly within the interval M=1.97±0.04 Mq as it has 
been reported for PSR J1614-2230. As we understand from the previous paragraphs 
it is possible to obtain more massive solutions by increasing rjv or decreasing i^jj. 
As both of these operations increase the transition density one eventually obtains 
massive NS which are purely hadronic. This can be either because the transition 
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Fig. 3. The overall result of the parameter study of the NJL model concerning the 
vector {rjv) and diquark (r;^)) channel coupling. The red band and all parameter 
pairs over it correspond to EoS which reproduce at least 1.97 ± 0.04 M0 (PSR 
J1614-2230), in the cyan region no stable hybrid configurations are found, while 
the grey hatched region corresponds to 'mostly' quark stars. A more detailed 
discussion is given in the text. Black filled circles correspond to parameterizations 
mentioned in the text. 

occurs at very large densities which are not realized in NS or, more relevant for 
our considerations, because the hybrid NS solutions become unstable. Examples 
for this situation are found in Fig. [5] for all r]v > 0.4. The cyan hatched region 
in Fig. [3] corresponds to all EoS parameterizations which result in unstable hybrid 
solutions. The opposite extreme scenario results from increasing r^o or decreasing 
r]v ■ Then, the transition density is lowered until eventually only a thin hadronic 
layer remains and the NS are basically pure quark star configurations. This scenario 
corresponds to the grey hatched region labeled 'Quark Stars'. In this parameter 
region the quark matter EoS has an early onset of the pressure due to a decrease 
of the dynamical quark mass before the first order phase transition which results 
in a low QM transition density. 

It is interesting to observe, that there is a band (red) of EoS parameters in the 
Vd-Vv plane with resulting maximum NS masses of 1.97 ±0.04 Mq. While the red 
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Fig. 4. Systematics of mass-radius and mass-central density curves for chosen pa- 
rameter sets with maximum mass equal to that of PSR J1614-2230. Large values 
of rjo and ryy result in a phase transition at lower densities and therefore a higher 
content of QM for massive NS. 

band in Fig. [3] denotes configurations with maximum masses corresponding to PSR 
J1614, it is worthwhile to ask for the actual quark content of these configurations. 
For the purpose of this article we will make a qualitative statement, only, and 
postpone quantitative analyses to a later publication. The curves labeled Monset 
in Fig. [3] indicate from which NS mass on the EoS results in hybrid NS solutions. 
Hence it is favorable to have a small value of Monset in order to obtain large QM 
cores. Therefore, the largest QM content for a high mass NS has to be expected in 
the upper right corner of Fig. [31 This is illustrated in Fig.Hl We point out, that the 
red band describes NS configurations where the maximum mass corresponds to the 
reported mass of PSR J1614, only. Of course, more massive NS are possible. As we 
are interested in massive NS with large quark cores Fig. [S] shows NS configurations 
from the upper right corner of Fig. |31 Large quark cores in our EoS parameteriza- 
tion are in general favored for values of rjy = 0.5 . . . 0.6. This is in vicinity of the 
value rjy = 0.5 as obtained after Fierz transformation of the one-gluon exchange 
interaction (details in Ref. [HI). It is noteworthy, that comparable small changes 
of the diquark coupling constant strongly affect the transition density. As an ex- 
ample. Fig. [S] illustrates, that with constant rjv = 0.6 a change of a few percent in 
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rjD (= 1-11 ■ • ■ 1.14) makes the difference whether a rather typical NS with a mass 
of, e.g., M= 1.4 M0 has QM content or not. 




Fig. 5. Mass-radius and mass-central density relation for a selection of EoS param- 
eterizations with mostly significant QM core (as found in the upper right corner of 
Fig. [31). 

3. 2. Implications for HIC 

In isospin-symmctric matter as found in HIC our main statements regarding 
the influence of r]v and iid on the stiffness of the EoS and the transition density do 
not change. Increasing rjy at constant rj^ increases the stiffness and the transition 
density, increasing jy^i at constant rjy reduces the stiffness and lowers the transition 
density. The transition densities in symmetric matter are plotted in Fig. [31 labeled 
as n^^^j. It is remarkable that along the red band (configurations with maximum 
masses corresponding to the mass of PSR J1614) the transition density in sym- 
metric matter has an almost constant value of n^^^^ « 4n5. This is a distinct 
difference to what we learned for the electrically neutral and /3-equilibrated EoS 
for NS matter, where the transition density decreased from the lower left to the 
upper right within the red band. Since the transition density in symmetric matter 
is lower only below the red band, a region where the maximum NS masses are not 
sufhcicntly large this model predicts n^^j^j w 4ns a-s a lower limit on the critical 
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density for the phase transition in symmetric matter. It is necessary to point out, 
that our previous study [3] gave a different result and predicted a lower transition 
density. In order to illustrate this, we plot hybrid EoS parameterizations in sym- 
metric matter corresponding to NS EoS with maximum masses of 1.97 ± 0.04 Mq 
(along the red band in Fig. ^ in Fig. IH] and our previous result in Fig. [T] While 
Fig. |6] suggests, that the phase transition occurs at densities too high to prevent 
the nuclear DBHF EoS from violating the flow constraint, our old results in Fig. [7] 
perform significantly better. This observation gives our study a surprising twist 
which wc will investigate in future work. The reason for this different behavior is 
not hard to comprehend. Both EoS start from a different parameterization already 
on the level of the scalar coupling constant. The larger dynamical quark mass 
resulting from the parameterization applied for the present study (see |35j for the 
actual values) causes a general shift of the QM onset to larger chemical potentials. 
Therefore, the phase transition from nuclear to QM simply follows this trend. From 
our perspective, this is the first time, that high density observables as elliptic flow 
and maximum NS masses can actually be applied to distinguish models with differ- 
ent parameterizations in the scalar channel. Of course, the usual caveats have to be 
made, lead by the most serious one: In order to draw final conclusions first it has to 
be made clear beyond any doubt, that NS actually have a QM core. At the current 
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Fig. 7. In previous work [Oj better agreement with the flow constraint has been 
achieved. Reason is a different parameterization which in particular resuhs in a 
smaller dynamical quark mass in vacuum. 



stage, this statement is neither proven nor disproven. If observational evidence in 
favor of the presence of QM in NS should ever arise, many more statements can 
be made, not only about the transition density in symmetric matter. Large NS 
masses require rather stiff QM EoS. As a consequence, the transition region in the 
density domain is not very large. This is illustrated Table [1] by a few examples for 
parameters within the red band in Fig. [3]) by calculating the difference An between 
QM and NM density at the phase transition. Since the quark content increases 
with rjv a clear signature for QM in NS would imply at least An<n5 Ri0.16 fm~'^. 



VD 


0.1 0.2 0.3 0.5 0.6 
0.92 0.98 1.02 1.10 1.15 


An [fm^"*] 


0.144 0.107 0.086 0.041 0.017 



Table 1. Density difference between QM and NM at the phase transition for a few 
parameterizations with M,nax G (1.97 ± O.O4)M0. 
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4. Conclusions 

We performed a systematic model analysis of an NJL-type EoS where we in- 
vestigated the dependence of NS properties on the diquark and vector coupling 
constants. It is possible to describe NS configurations with a significant amount of 
QM and maximum NS masses which are in agreement with the reported high mass 
of PSR J1614-2230. The amount of QM in a NS will increase with increasing vector 
coupling if the diquark coupling is adjusted accordingly. Therefore, both channels 
are important for the understanding of NS phenomenology. Due to a different 
choice of the scalar coupling and the resulting larger dynamical quark mass the 
present model parameterization does not resolve the conflict of the nuclear DBHF 
EoS with the flow constraint. This result is of great importance since it implies 
that without a sound understanding of the chiral phase transition in medium oth- 
erwise identical models can come to quantitatively significantly different results. 
In particular this raises doubts about the predictive power of any EoS which does 
not account for the mechanism of chiral symmetry breaking. Of course, we do not 
intend to imply, that NJL-type models will provide the ultimate tool to describe 
the QCD-phase transition qualitatively and quantitatively correct. There are many 
open questions, which will need to be addressed in the future, most importantly 
the mechanism of confinement / deconfinement in medium. Significant improvement 
of any existing model is required in order to understand this phenomenon and to 
finally get rid of the necessity to construct thermodynamically motivated Maxwell- 
and Gibbs-phase transitions, which are not suited to provide deeper insights into 
the microphysical mechanisms governing QCD phase transformations. First steps 
in this direction have been performed within a generalized N JL model |43j . Further 
progress was made in a microscopic description of the baryon dissociation in dense 
matter, based again on a NJL model approach |44| where in particular the interplay 
between chiral symmetry restoration and diquark condensation transitions at high 
densities for the spectral function of nucleons has been investigated. Further steps 
will be taken towards the development of an EoS which describes the Mott disso- 
ciation of baryons into their quark constituents in a consistent way which accounts 
for nucleonic bound and scattering states simultaneously. An EoS of this quality 
has been developed within a generalized Beth-Uhlenbeck approach for strongly in- 
teracting matter in order to describe the Mott dissociation of deuterons in nuclear 
matter |45] (see also recent work on cluster formation in low-density nuclear matter 
[151 ITT] ) and the Mott effect for mesons in quark matter |48l . We are convinced 
that studies of the Mott mechanism for the dissociation of baryons in dense mat- 
ter which generalize early nonrelativistic approaches [50l [51] to a field theoretical 
formulation will illuminate the role of this mechanism in its interplay with chiral 
symmetry restoration and diquark condensation for the EoS of dense matter and 
its applications in HIC and in astrophysical environments. 
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